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Abstract 

We investigate the 2D ferromagnetic Kondo lattice model for manganites with classical corespins at Hund's rule 
coupling Jh = 6, with antiferromagnetic superexchange 0.03 < J' < 0.05. We employ canonical and grand canonical 
unbiased Monte Carlo simulations and find paramagnetism, weak ferromagnetism and the Flux phase, depending 
on doping and on J' . The observed single particle spectrum in the flux phase differs from the idealized infinite 
lattice case, but agrees well with an idealized finite lattice case with thermal fluctuations. 
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The ferromagnetic Kondo lattice model is believed 
to describe the electronic degrees of freedom of com- 
pounds such as manganites. In these materials the spins 
of the itinerant electrons have a strong ferromagnetic 
coupling to localized corespins Si, leading to an Effec- 
tive Spinless Fermion (ESF) Hamiltonian [I] 



H = 
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where c\ creates an electron at site i with its spin par- 
allel to Si, and Jh = 6 throughout this work. The 
hopping of these spinless electrons is modulated by 
the factor t] ■ = cos(l?ij/2) e 1 ** 3 , which depends on 
the relative angles #jj of the two neighboring corespins 
and carries a complex phase «/>y . On a two-dimensional 
lattice, this complex phase leads to the so call "Flux 
phase" around half filling of the lower Kondo band 
(doping x = 50%) [2]. While the system is ferromag- 
netic at half filling for J' < 0.025, the Flux phase has 
been shown [2] to be the ground state at half filling for 
larger values of J' . 
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We performed canonical and grand canonical unbi- 
ased Monte Carlo simulations [3] . The magnetic phases 
can be observed in the corespin structure factor which 
is depicted in Fig. I for J' = 0.05. In the polaronic 
regime with high compressibility at x « 18% holes [3], 
only an AFM signal is visible. At x w 24%, the com- 
pressibility is suddenly much reduced and the system 
becomes homogeneous but a small AFM peak persists. 
At x fa 28% no magnetic structure at all is visible, 
The flux phase appears at x « 44% and x = 50% with 
signals at (tt,0) and (0, 7r). The only FM signal exists, 
weakly, at even larger doping x ~ 62%. 

Fig. 2 shows the normalized FM corespin moment 
S(q = 0) = | Ytii Si\/ L as a function of doping, where L 
is the number of lattice sites. It is expected to vanish in 
the Flux phase. One sees that the doping range of the 
Flux phase becomes smaller at smaller J' . The lattice 
is weakly ferromagnetic at slightly lower doping x w 
40%. All our data are compatible with the existence of 
a flux phase in the region where S(q — 0) ~ 0. 

The corespin structure in an idealized flux phase is 
depicted in Fig. 3(a). While all corespins are drawn par- 
allel to the xy-plaxie for better visibility, a global rota- 
tion in spin space leads to equivalent configurations. In 
such corespin configurations, \t\ A = \j\p2 is constant. 
However, the hopping carries a complex phase tpij. Di- 
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Fig. 1. Corespin structure factor for Jh — 6, & — 50, J — 0.05 
and various doping levels. Simulations were done on a 12 X 14 
lattice, i. c. for L — 168 lattice sites. 
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Fig. 2. FM corespin Moment S(q — 0) as a function of doping 
for various values of J', remaining parameters as Fig. 1. 

agonalizing the ESF Hamiltonian for the four site unit 
cell of the Flux phase leads to the dispersion relation 



e(k x ,k y ) = — -i ± \/2(cos 2 k x + cos 2 k y ) 



(2) 



The one-particle DOS of these bands for the infinite 
lattice with the corespins aligned in the perfect Flux 
phase is depicted in Fig. 3(b); it features a pseudogap 
at half filling. The dispersion and the DOS of the full 
Kondo model and for Jh — ► oo can be found in [2] . 

The spectral density of the Flux phase is subject 
to strong finite size effects, which are caused by the 
fact that only a very small part of the Brillouin zone 
contributes to the few states around the pseudogap at 
(•7r/2, 7r/2). In order to ascertain that we did indeed ob- 
serve the Flux phase, we therefore compared our MC 
results to a model system. We started from a perfect 
Flux phase configuration of the corespins on a finite 
lattice, added random perturbations, and diagonalized 
the resulting ESF Hamiltonian. The amount of fluctua- 
tions was chosen to fit the width of the peaks in the full 
MC results. Fig. 4 compares the average spectral den- 
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Fig. 3. Perfect Flux phase: (a): Schematic plot of the corespin 
configuration. The dashed square delimits the unit cell. All 
spins arc drawn parallel to the icy-plane, (b): Electronic DOS 
for the ESF Hamiltonian with Jh — 6 on an infinite lattice. 
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Fig. 4. Spectral density for the Flux phase at half filling for 
j' = 0.05; P = 50 (left) and (3 - 80 (right) on a 14 X 12-lattice. 
Dashed lines: MC results for the ESF model. Solid lines: ideal 
Flux phase with random deviations. 

sity and the DOS thus obtained to MC simulations for 
the ESF model exhibiting almost perfect agreement, 
while there is a notable difference between the DOS in 
Fig. 4 and the idealized case in Fig. 3(b), The system 
can thus indeed be well described by thermal fluctua- 
tions about the Flux phase on a finite lattice. 
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